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The external-field QCD Sum Rules method is used to evaluate the coupling constant of the light 
, isoscalar-scalar meson ("cr" or e) to the nucleon. The contributions that come from the excited 

nucleon states and the response of the continuum threshold to the external field are calculated. The 
obtained value of the coupling constant is compatible with the large value required in one-boson 
£\J , exchange potential models of the two-nucleon interaction. 
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l I I. INTRODUCTION 

> : 

. The values of the meson-baryon couplingconstants are of particular interest in understanding the nucleon-nuclcon 
(NN) 0, Q and hyperon-nucleon (YN) 0, Q interactions in terms of e.g. one-boson exchange (OBE) models. The 
scalar mesons play a significant role in such phcnomcnological potential models. The structure and even the status 
of the scalar mesons have, however, always been controversial . In early OBE models for the NN interaction 
the exchange of an isoscalar-scalar "cr" meson with a mass of about 500 MeV was needed to obtain enough medium- 
i range attraction and a sufficiently strong spin-orbit force. R was only later understood that the exchange of a broad 
isoscalar-scalar meson, the e(760), simulates the exchange of such a low-mass "cr" 0. Thee(760) is difficult to detect 
because it is broad and hidden under the strong signal from the p°(770). There are strong arguments from chiral 
,_L| ' symmetry for the existence of such a light isoscalar-scalar meson approximately degenerate with the p meson Q . 

In the quark model, the simplest assumption for the structure of the scalar mesons is the 3 Po qq states. In this case, 
the scalar mesons might form a complete nonet of dressed qq states, resulting from e.g. the coupling of the P-wave 
qq states to meson- meson channels [9|- Explicitly, the unitary singlet and octet states, denoted respectively by E\ and 



o 
m 
O 



=5 



- 1—1 

X 



£8, read 



£\ = (uu + dd + ss) I V3 , 



e 8 = (uu + dd-2ss)/VE. (1) 
The physical states are mixtures of the pure SU (3)-fiavor states, and are written as 

e = cos 9 s Si + sin S e 8 , 

fo = - sin 6 S si + cos 9 S e 8 . (2) 
For ideal mixing holds that tan# s = 1/^/2 or 9 S ~ 35.3°, and thus one would identify 

e(760) = (uu + dd)/V2, 

/o(980) = -aa. (3) 
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The isotriplet member of the octet is aj'°(980), where 

ag(980) = (uu-dd)/V2. (4) 

An alternative and arguably more natural explanation for the masses and decay properties of the lightest scalar 
mesons is to regard these as cryptoexotic q 2 q 2 states [ic| . In the MIT bag model, the scalar qq states are predicted 
around 1250 MeV, while the attractive color-magnetic force results in a low- lying nonet of scalar q 2 q 2 mesons Jiol lllj . 
This nonet contains a nearly degenerate set of 7 = and 1=1 states, which are identified as the / (980) and a ' (980) 
at the KK threshold, where 

a°(980) = {sdsd- susu)/V2, 

/ (980) = (sdsd + susu)/V2, (5) 

with the ideal-mixing ang le tan0 s = -a/2 or 9 S ~ -54.8° in this case. The light isoscalar member of the nonet is 

e(760) = udud. (6) 

The nonet is completed by the strange member k(880), which like the e(760) is difficult to detect because it is hidden 
under the strong signal from the if* (892) 0,0. In keeping with other recent works 0, El 0] we will use in this 
paper the nomenclature (a^'° , fa, cr, n) for the scalar-meson nonet, where one should identify a = e(760), but we will 
not rely on a particular theoretical prejudice about the quark structure of the light scalar mesons. 

One way to make progress with the scalar mesons is to study their role in the various two-baryon reactions (NN, 
YN, YY). Our aim in this paper is to calculate the nucleon-a coupling constant gNNa by using the QCD Sum Rules 
(QCDSR) method |l5|. QCDSR links the hadronic degrees of freedom with the underlying QCD parameters, and 
serves as a powerful tool to extract qualitative and quantitative information about hadron properties [Trl Il7| . In this 
framework, one starts with a correlation function that is constructed in terms of hadron interpolating fields. On the 
theoretical side, the correlation function is calculated using the Operator Product Expansion (OPE) in the Euclidian 
region. This correlation function is matched with an Ansatz that is introduced in terms of hadronic degrees of freedom 
on the phenomenological side. The matching provides a determination of hadronic parameters like baryon masses, 
magnetic moments, coupling constants of hadrons, and so on. 

The QCDSR method has been extensively used to investigate meson-baryon coupling constants. One usually starts 
with the vacuum-to-vacuum matrix element of the correlation function that is constructed with the interpolating 
fields of two baryons and one meson. However, this three-point function method has as a major drawback that at 
low momentum transfer the OPE fails. Moreover, when the momentum of the meson is large, the latter is plagued 
by problems with higher resonance contamination 18]. A method that can be used at low momentum transfer is the 
external-field method [l^. There are two formulations that can be used to construct the external-field sum rules: The 
first one is to start with a vacuum-to-vacuum transition matrix element of the nucleon interpolating fields. In this 
approach, no vacuum-to-meson matrix elements occur, but one has to know the response of the various condensates 
in the vacuum to the external field, which can be described by a susceptibility \- This method has been used to 
determine the magnetic moments of baryons |l9j, l20l l2lL |22j , the nucleon axial coupling constant [22t |23| , the nucleon 
sigma term |24j . and baryon isospin mass splittings |25| . In the second approach, one starts with a vacuum-to-meson 
transition matrix element of the nucleon interpolating fields, where some other transition matrix elements should be 
evaluated 16]. (This is also the starting point of the light-cone QCDSR method.) In 26], pion-nucleon coupling 
constant was calculated in the soft meson limit using this approach. Later it was pointed out that the sum rule 
for pion-nucleon coupling in the soft-meson limit can be reduced to the sum rule for the nucleon mass by a chiral 
rotation so the coupling was calculated again with a finite meson momentum |27j . These calculations were improved 
considering the coupling schemes at different Dirac structures and beyond the chiral limit contributions [23, |2j|, |3(| • 
This coupling constant has also been calculated using the vacuum-to- vacuum method [3ll l32l] , and it was pointed out 
that the sum rule that is constructed for the coupling is independent and it is not reduced to the nucleon mass sum 
rule by a chiral rotation. 

In this paper, we calculate the nucleon-er coupling constant g^Na by using the external-field QCDSR method. 
We evaluate the vacuum-to-vacuum transition matrix element of the two-nucleon interpolating fields in an external 
isoscalar-scalar field, and construct two sum rules, one of which leads to a stable result with respect to variations 
in the Borel mass. We also compute the contributions that come from the excited nucleon states and the response 
of the continuum threshold to the external field. Previously, the strong and weak (parity-violating) pion-nucleon 
coupling constants [3ll 13^ ] and the coupling constants of the vector mesons p(770) and w(782) to the nucleon [34| 
were calculated by using this method. 

We will compare our result for the coupling constant with the value from a successful OBE model of the NN 



ig con 

interaction, the Nijmegen soft-core potential [1 0, which was originally derived from Regge-pole theory. (The 
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coupling constants of this OBE model were analyzed from the point of view of the large- N c expansion of QCD in 
Ref. |35j.) It is then important to realize that in NN potential models the coupling constants of the heavy mesons to 
the nucleon are determined by the ("non-peripheral") S-, P-, and I?- waves. Therefore, the fits to the scattering data 
are sensitive to e.g. the volume integral of the OBE potentials, which is proportional to the coupling at t — j3|| 
(< = — p 2 , where p is the four-momentum of the meson). The coupling constants obtained from the external-field 
QCDSR method are also defined at t = 0, and therefore the comparison to the OBE model is appropriate. 

Our paper is organized as follows: In Section[H]we present the formulation of QCDSR with an external isoscalar- 
scalar field and construct the relevant sum rules. We give the numerical analysis of the sum rules and discuss the 
results in Section ITTT1 Finally, we arrive at our conclusions in Section llVl 



II. NUCLEON SUM RULES IN AN EXTERNAL SCALAR FIELD 



In the external-field QCDSR method one starts with the correlation function of the nucleon interpolating fields in 
the presence of an external constant isoscalar-scalar field er, defined by 



11° (q) 



d 4 x e*' (o\T[ m (x)fj N (0)} 



where rj^ is the Ioffe nucleon interpolating field |33 



(7) 



(8) 



Here a, b, c denote the color indices, and T and C denote transposition and charge conjugation, respectively. The 
external scalar field contributes to the correlation function in Eq. in two ways: First, it directly couples the quark 
field in the nucleon current and second, it modifies the condensates by polarizing the QCD vacuum. In the presence 
of an external scalar field there are no correlators that break Lorentz invariance, like (qcr^q) which appears in the 
case of an external electromagnetic field F^ v . However, the correlators already existing in the vacuum are modified 
by the external field, viz. 



(qq) a = (qq) + g a q X^{m) , 
(g c qa- ■ Gq) a = (g c q<T ■ Gq) + g q XG^{9 c qa- ■ Gq) , 



(9) 



where g^ is the quark-cr coupling constant and, x an d Xg are the susceptibilities corresponding to quark and mixed 
quark-gluon condensates, respectively. We have assumed that the responses of the up and the down quarks to the 
external (isoscalar) field are the same. 

In the Euclidian region, the OPE of the product of two interpolating fields can be written as 



(10) 



where C°(q) are the Wilson coefficients and O n are the local operators in terms of quark and gluon fields. At the 
quark level, we have 



T[ m (x)fj N (0)] 



2i ^f'MT [ {S*(x)i u C[S?\x)] T Ci lt }>r 6 'fS?(x)>fTB- 



(11) 



In order to calculate the Wilson coefficients, we need the quark propagator in the presence of the external sigma field. 
In coordinate space the full quark propagator takes the form 



S q (x)=S^(x) + S^(x), 



(12) 



where 



i S^ ab 



(0\T[q a (x)q b (0)\0) o 



b yc 



2tt 2 x 4 ' 



G 1 ' 



1 



32tt 2 2 ^x 1 



(a^x + xa^) 



5 ab 

~12 



S ab x 2 

192 



(g c qa- ■ Gq) , 



(13) 
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and 



i S q a)ah 



(0\T[q a (x)<?(0)\0} a 

gab 1 



9 a O 



47r 2 a 
i S ab 



32tt 



2 . S ab (g 2 G 2 / 2 



S ab x 2 

192 



12 



XG(g c q<r ■ Gq) 



i S ab x 2 
2 7 x 3 2 



2 9 x 3tt 2 



x 2 ln(-a; 2 ) 



(14) 



In these expressions, G^ v is the gluon field tensor and g 2 = 4ira s is the quark-gluon coupling constant squared. We 
do not include terms that are proportional to the quark masses, since these terms give negligible contributions to the 
final result. 

Using the quark propagator in Eq. l|12f) . one can compute the correlation function II <T (g). The diagrams that we 
use to calculate the Wilson coefficients of W (q) are shown in Fig. ^ Lorentz covariance and parity imply that IF (q) 
takes the form 



(15) 



where q — <z Al 7 AI . Here IlJ and IIq represent the invariant functions in the vicinity of the external field, which can be 
used to construct the mass sum rules for the nucleon, and II* and denote the invariant functions in the presence 
of the external field. Using these invariant functions, one can derive the sum rules at the structures 1 and q. H% and 
III are evaluated as follows: 



n»Gz) 



'« (27T) 



-q ) - X ^ % + —a q - (x + XG) ^4 a q 



(16) 



and 



(27T) 



9 4 l / 2\ 10 2 2 2 TTlg 2 \ 

y ln(-g ) - -^a q - \a q q m(-g ) + XGya, m(-g ) 



ln(-g 2 



X 



24q 



2 u. q 



(17) 



(2^) 2 (^), b = ( 5c 2 G 2 >, and m 2 = (g c q* ■ Gq)/(qq). 



where we have defined a q ; 

We now turn to the calculation of the hadronic side. We saturate the correlator in Eq. (Q) with nucleon states and 
write 



n- (?) = <*^<^ w) f'»l°> 



N 



M 2 ' 



(18) 



where M/v is the mass of the nucleon. The matrix element of the current t]n between the vacuum and the nucleon 
state is defined as 



\vn\N) = X N v , 



(19) 



where Ajv is the overlap amplitude and v is the Dirac spinor for the nucleon, normalized as vv = 2M/\r. Inserting 
Eq. lO into Eq. (HE) and defining gNNa via the interaction Lagrangian density 



-gNNa vva 



(20) 



we obtain for the hadronic part 



,2 q + M N 



gNNo 



q + M N 



q*-M" N q* 

In addition, there are contributions coming from the excitations to higher nucleon states which are written as 

. . q + M N q + M N 

AjvAat. — ^-py- gNN» 



(21) 



M 2 N 



2 -M 2 N » 



(22) 
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as well as contributions coming from the intermediate states due to a-N scattering, i.e. the continuum contributions. 
The term that corresponds to the excitations to higher nucleon states also has a pole at the nucleon mass, but a single 
pole instead of a double one like in Eq. (|21(1 . This single-pole term is not "damped" after the Borel transformation 
and should be included in the calculations. 

Finally, there is another contribution that comes from the response of the continuum to the external field, given by 



-As 6(g) w 

- — d(s - s )ds . 

s — q z 



(23) 



where so is the continuum threshold, Asq is the response of the continuum threshold to the external field, and b(s) 
is a function that is calculated from the OPE. When Aso is large, this term should also be included in the hadronic 
part EH. 

The QCD sum rules are obtained by matching the OPE side with the hadronic side and applying the Borel 
transformation. The resulting sum rules are: 



-M*a q E + ^ X M*alL^ 



^ M 2 a q _ (x + Xg) «£ a 2 L - 2/ 27 



~ A N 



M 



N 



■ gNNo 



_ Ml! 

9 a q ^ 2 fl J 



As q Q M 2 L~ 4 / 9 e( M "- s °V M 



(24) 



and 



20 



2 M 8 E 2 i" 4/9 + — M 2 a 2 L 4/9 + 2 x a q M 6 E 1 



XG m a q 



o a q M h,q h 



— - M Eq L 



-4/9 



-{2M 



12 

M 2 ) — 



■ M 2 a c 



gNNcr 



where M is the Borel mass and we have defined \ 2 N 
factors 



- M 2 
f-Bi — 

32tt 4 X 2 n . 



(25) 



The continuum contributions are included by the 



E Q = l-e- s ^ M2 



E x = l-e^o/M 



E 2 = l-e 



i 1 ' M 2 J 
V M 2 2A/ 4 J 



(26) 



where Sq are the continuum thresholds with i = l,q. In the sum rules above, we have included the single-pole 
contributions with the factors Bi. The third terms on the right-hand side (RHS) of Eqs. I|24|) and l|25|) denote the 
contributions that are explained in Eq. i|23[l . These terms are suppressed by the factor e - ( s o _M «)/ M2 as compared to 
the single-pole terms. We have incorporated the effects of the anomalous dimensions of the various operators through 
the factor L 



HM 2 /A 2 QCD 



)/Hn 2 /K 2 QCD ). 



III. ANALYSIS OF THE SUM RULES AND DISCUSSION 

In this Section we analyze the sum rules derived in the previous Section in order to determine the value of gNNa ■ 
We observe that the sum rule in Eq. \2ty is more stable than the other sum rule in Eq. (|25ll . so we use only this sum 
rule for the numerical analysis. Such a comparison and conclusion have been made about these sum rules also in some 
earlier works Ell- 
in order to calculate gNNa- , we need to know the values of the scalar susceptibilities \ and xg- The value of the 
susceptibility x can be calculated by using the two-point function |2J| 



T(p 2 )=i / d 4 xe ip - x (0 T[u(x)u(x) + d(x)d(x),u(0)u(0) + d(0)d(0)] 0) , 



(27) 



via the relation 



X (qq) = =2X0). 



(28) 
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The two-point function in Eq. I|27[l at p 2 = has been calculated in chiral perturbation theory (3i| with the result 

(W> fh , 7 , 11' 



16^ 2 / 4 



(^ + 5«"-t)' (29) 



where f % = 93 MeV is the pion decay constant and i\ and £2 are low-energy constants appearing in the effective 
chiral Lagrangian. The values of these low-energy constants have been estimated previously in various works (see e.g. 
Ref. ^(J f° r a review). A recent analysis of it-it scattering gives i\ — —1.9 ± 0.2 and £2 — 5.25 ± 0.04 [4jj, which 
is consistent with earlier determinations, but with smaller uncertainties. Using these values of £\ and £2 and taking 
the quark condensate a q = 0.51 ± 0.03 GeV 3 , we find \ = —10 ± 1 GeV -1 . The value of the susceptibility \G is less 
certain. Therefore, we allow xg to vary in a wider range. We also adopt b — 4.7 x 10 -1 GeV 4 , = 2.1 GeV 6 , 
and m 2 = 0.8 GeV 2 [H El . We take m n = 0.94 GeV, the renormalization scale /1 = 0.5 GeV, and the QCD scale 
parameter Aqcd = 0.1 GeV. It is relevant to point out that the choice of the two-point function in Eq. 1(2 7|) does not 
imply a theoretical prejudice about the structure of the scalar mesons. What is calculated is just the susceptibility 
pertaining to the response of the quark condensates (qq) to the scalar qq field, as shown in Eq. 
To proceed to the numerical analysis, we arrange the RHS of Eq. 1(24(1 in the form 

f(M 2 ) = A q + B q M 2 + C q M 2 L-^e^- s ^l M °- , (30) 

and fit the left-hand side (LHS) to f(M 2 ). Here we have defined 

a — \2 M n 

A q = —*N — —9NNa, 

9q 



B q 



9, 



C q = ^Asl (31) 



In Fig. |3 we present the Borel mass dependence of the LHS and the RHS of Eq. I(24() for Sq = 2.3 and xg = X = — 10 
GeV -1 . We choose the Borel window 0.8 GeV 2 < M 2 < 1.4 GeV 2 which is commonly identified as the fiducial region 
for the nucleon mass sum rules 

EH- 

It is seen that the LHS curve (solid) overlies the RHS curve (dashed) . In order to 
estimate the contributions that come from the excited nucleon states and the response of the continuum threshold, we 
plot each term on the RHS individually. We observe that the single-pole terms (dotted) give very small contribution, 
but the response of the continuum threshold (dot-dashed) is quite sizable. Nevertheless, the summation of these 
curves with the line of the double-pole term (small-dashed) gives a stable sum rule. 

In Fig. |2| we plot the Borel mass dependence of the four terms on the LHS of Eq. (|24|l separately, together with 
their summation for Sq = 2.3 GeV 2 and xg = X = ~ 10 GeV -1 . This helps us to compare the contributions of different 
operators on the OPE side. Here 0\ denotes the first term, 2 denotes the second term, and so on. We observe that 
0i and 03 are small, O4 is sizable, and 02 is large. The term O4 contributes with different sign with respect to 0\ 
and 03, and so tends to cancel the latter. Therefore gNNa is mainly determined by 02 on the LHS. 

In order to see the sensitivity of the coupling constant on the continuum threshold and the susceptibility %, we plot 
in Fig.^the dependence of gNNa/g q on X f° r three different values Sq = 2.0, 2.3, and 2.5 GeV 2 , and taking x = Xg- 
One sees that gNNa- changes by approximately 8% in the considered region of the susceptibility x- The value of gNNa- 
is not very sensitive to a change in Sq, which gives an uncertainty of approximately 6% to the final value. Taking into 
account the uncertainty in x, Sq, and a q , the predicted value for gNNa/ g q of the sum rule in Eq. <|24|l reads 

9NN*/9 q =3.9 ±1.0. (32) 

In a similar way, one can calculate the other two terms on the RHS of Eq. (|24|l as: 

B q = -0.2 ± 1.2 GeV 5 , 

C q = -7.9 ± 2.9 GeV 5 . (33) 

As noted above, the value of the susceptibility xg is less certain than the value of %. If we let xg change in a wider 
range, say 6 GeV -1 < — xg < 14 GeV -1 , this brings an additional 15% uncertainty to the value quoted in Eq. Ipffil, 
The ratio in Eq. 1)32(1 is in agreement with the naive quark model, which gives gNNa/ g q = 3 based on counting the 
u- and the d-quarks in the nucleon. (Ideal mixing in the scalar sector is assumed above, that is, the sigma meson is 
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taken without a strange-quark content.) Another estimate can be made from the ratio of pion-nucleon to pion-quark 
coupling constant, gNNn/gq- Since the a meson is the chiral partner of the pion |8|, one expects 

9NNa/9q = 9NN-n/gl ■ (34) 

Using the Goldberger-Treiman relation for both the pion-nucleon and the pion-constituent quark couplings, 

A M N 
9nntt = 9n —f~ > 

Jit 

9^ = 9^, (35) 

Jit 

where m q is the mass of the constituent quark, g N and g q are the nucleon and the quark axial couplings, respectively, 
one obtains |4^ | 

gNN-rr _ 5 M N 
Sj ~ 3 m, ' 

With a constituent-quark mass of 340 MeV H2, Eq. JSBJl yields gNNn/gq = 4.6. Using Eq. (P3J we find that this 
agrees nicely with the QCDSR result in Eq. I|32ll . 

To determine gNNa , one next has to assume some value for the quark-er coupling constant g a q . Adopting the value 
g° = 3.7 as estimated from the sigma model I43j . we obtain 

g NNa = 14.4 ± 3.7 . (37) 

The coupling constant in Eq. (|3TJl is defined at t = 0, i.e. gNNa = gNNa(t = 0). As stressed above, also in NN 
potential models the heavy- meson coupling constants are determined at t — 0. The (large) value of gNNa obtained 
in Eq. (|o7[) is in agreement with the value gNNa — 16.9 from the Nijmegen soft-core NN potential model pj, obtained 
from a fit to the NN scattering data. 



IV. CONCLUSION 



We have calculated the coupling constant gNNa of the isoscalar-scalar meson, which plays a significant role in OBE 
models of the NN and YN interactions, to the nucleon, using the external-field QCDSR method. Our main result is 
the ratio gNNal g a q in Eq. (|32|l which is determined purely from QCDSR. The value of gNNa is dependent on g a q , the 
value of which we use as estimated in the sigma model. The obtained value of gNNa is in agreement with the large 
value found in OBE models. We have also computed the contributions that come from the excited nucleon states and 
the response of the continuum threshold to the external field. We observe that while the single-pole contributions are 
small, the response of the continuum threshold is sizable. We plan to extend the external-field QCDSR method to 
the hyperons and the complete scalar-meson nonet, in order to address the S't r (3)-flavor structure of the scalar-meson 
coupling constants to the baryon octet [Sj . 
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FIG. 1: The diagrams that were used to calculate the Wilson coefficients of the correlation functions 11^ and n*. The solid, 
wavy, and the dashed lines represent the quark, gluon, and the external scalar field, respectively. 



10 





-5 

^ -10 
-15 





1 


1 






- 









C q M 2 e (M N- 


LHS 

RHS - 

A, 

B„M 2 - 

i 


i 



0.8 1 1.2 1.4 

M 2 (GeV 2 ) 

FIG. 2: (Color online) The Borel mass dependence of LHS and the fitted RHS of Eq. ^ for = 2.3 GeV 2 and X a = X = -10 
GeV -1 . We also present the terms on the RHS individually. Note that the LHS curve (solid) overlies the RHS curve (dashed). 
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FIG. 3: (Color online) The four terms on the LHS of Eq. 1241 individually, together with the summation of them for Sg = 2.3 
GeV 2 and xg = X — — 10 GeV -1 . Here Ox denotes the first term, O2 denotes the second term, and so on. 
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FIG. 4: (Color online) The dependence of gjvjvtx/ff? on the susceptibility x f° r three different values of Sq = 2.0, 2.3, and 2.5 
GeV 2 ; here we take \ = XG- 



